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>infrasonic  signal-to-noise  ratio  will  be  greater  at  the  output  of  a  vertical 
seismograph  than  at  the  output  of  a  microbarograph  in  a  frequency  range 
extending  roughly  from  0.005  Hz  to  0.05  Hz  during  intervals  of  atmospheric 
turbulence.  /^Further  research  activities  will  include  documentation  of  the 
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turbulence  intensity;  a  renewed  investigation  into  the  relationship  between 
the  yield  of  latmospheric  nuclear  explosions  and  the  spectrum  of  the 
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seismograph  system. 
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1.0  BACKGROUND 


There  Is  a  need  to  improve  upon  our  current  capability  to  detect  and  identify 
atmospheric  nuclear  tests.  In  particular,  it  appears  that  circumstances  can 
arise  where  considerable  uncertainty  may  surround  the  origin  and  source 
location  of  anomalous  signals  detected  by  the  existing  atmospheric  sur¬ 
veillance  systems.  It  is  generally  agreed  that  the  uncertainty  could  be 
significantly  reduced  with  supplemental  information  provided  by  an  infrasonic 
monitoring  system.  The  purpose  of  this  research  program  is  to  investigate  a 
new  method  of  acquiring  the  desired  infrasonic  information. 

Atmospheric  nuclear  explosions  produce  pressure  pulses  which  are  potentially 
detectable  at  long  ranges  from  the  source.  The  principal  energy  of  these 
pulses  has  been  found  to  lie  in  the  0.001  -  1  Hz  bandwidths.  In  the  past, 
the  capability  to  detect  the  infrasonic  signals  from  atmospheric  nuclear 
tests  was  vested  in  a  network  of  microbarograph  arrays.  These  arrays,  which 
now  have  been  largely  abandoned,  consisted  of  four  or  more  sensors  which  had 
the  intrinsic  sensitivity  to  resolve  atmospheric  pressure  changes  as  small  as 
a  few  hundredths  of  a  microbar  in  the  bandwidth  of  Interest  for  nuclear  test 
monitoring.  Unfortunately,  this  sensitivity  could  rarely  be  used  to  prac¬ 
tical  advantage  to  detect  infrasonic  signals  because  of  the  relatively  large 
atmospheric  pressure  fluctuations  intermittently  created  by  the  local  winds. 

In  order  to  combat  the  wind  noise  problem,  a  "pipe  array”  (cf.  Daniels,  1959) 
was  attached  to  the  inlet  of  the  microbarograph.  This  device  is  a  passive 
wavelength  filter  which  consists  of  a  hollow  pipe,  less  than  3  cm  in 
diameter,  with  maximum  linear  dimensions  on  the  order  of  several  hundred 
meters  or  more,  and  pressure  inlet  ports  placed  at  2-4  meter  intervals. 

Given  the  assumption  that  the  wind  noise  is  uncorrelated  between  inlet  ports 
while  the  Infrasonic  signal  remains  correlated,  it  is  then  possible  to  obtain 
an  improvement  in  infrasonic  SNR  approaching  where  N  is  the  number  of 
inlets.  Experiments  demonstrated  that  while  such  Improvements  were  indeed 
possible  at  relatively  high  frequencies  (say  greater  than  .2  -  .3  Hz),  the 
state  of  organization  of  the  wind-generated  noise  generally  precluded 
achievement  of  significant  infrasonic  SNR  Improvement  at  lower  frequencies 
(MacDonald  et.  al . ,  1971).  Thus,  the  "pipe  array"  is  of  limited  value  if  one 
is  interested  in  enhancing  Infrasonic  SNR's  over  the  entire  frequency  range 
of  interest.  Another  approach  vrtilch  was  under  serious  consideration  in  the 
early  1970's  was  to  increase  the  number  of  microbarographs  in  the  arrays  by 
at  least  a  factor  of  3  and  to  rely  on  multichannel  processing  techniques  to 
improve  the  Infrasonic  SNR  (cf.  Teledyne  Geotech,  1971).  While  implementation 
of  this  approach  is  certainly  well  within  the  current  state  of  the  art,  the 
probable  cost  of  developing  and  operating  an  advanced  infrasonic  network 
leads  to  the  consideration  of  alternative  methods. 

During  the  early  1970 's,  the  AF0SR  sponsored  basic  studies  of  the  rela¬ 
tionship  between  local  atmospheric  pressure  changes  and  low-frequency  seismic 
noise  (Sorrells,  1971;  Sorrells  and  Goforth,  1973;  Sorrells  and  Douze  1974; 
Savlno  et  al,  1972;  Savlno  and  Rynn,  1972).  An  Important  result  of  these 
investigations  was  the  prediction  and  experimental  confirmation  that  the 
quasi-static  earth  movements  triggered  by  the  passage  of  infrasonic  waves 
could  be  detected  at  the  outputs  of  a  sensitive  long-period  seismograph 
system.  Equally  important  was  the  prediction  and  observation  that  the  earth 
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acts  as  a  passive  wavelength  filter  with  respect  to  the  atmospheric  pressure 
fields,  selectively  attenuating  the  shorter  wavelength  components  to  a 
greater  degree  than  the  longer  wavelength  components.  The  properties  of 
this  filter  are  controlled  by  the  local  distribution  of  elastic  constants 
and  the  depth  of  observation  and  are  virtually  independent  of  the  state  of 
organization  of  the  input  field.  Since,  for  a  given  frequency,  the  convec¬ 
tive  wavelength  of  wind-generated  pressure  noise  may  be  an  order  of  magni¬ 
tude  shorter  than  that  of  an  infrasonlc  signal,  relatively  large  SNR  gains 
are  theoretically  possible  regardless  of  the  state  of  organization  of  the 
signal  and  noise  field.  This  point  is  illustrated  by  the  results  shown  in 
figure  1.  This  is  the  theoretical  infrasonlc  SNR  gain,  predicted  for  obser¬ 
vations  made  at  a  depth  of  100  meters  in  a  homogenous,  Isotropic,  perfectly 
elastic  half  space  with  equivalent  Lame'  constants.  The  results  shown  are 
appropriate  for  a  vertical  seismograph  and  assume  that  the  propagation  speed 
of  the  signal  is  330  m/sec  while  the  convection  velocity  of  the  wind¬ 
generated  noise  is  10  m/sec.  Observe  that  the  predicted  SNR  gain  is  greater 
than  30  dB  throughout  the  entire  frequency  range  of  interest.  Conceptually 
then,  the  earth-seismograph  system  may  be  thought  of  as  a  functional  equiva¬ 
lent  to  the  pipe  array-microbarograph  system  with  the  added  advantage  that 
the  response  of  the  former  system  is  virtually  independent  of  the  state  of 
organization  of  the  signal  and  noise  fields.  Thus,  a  high-quality,  long- 
period  seismograph  system,  coupled  with  a  single  microbarograph  to  aid  in  the 
discrimination  between  seismic  and  atmospherically  Induced  earth  movements, 
could  conceivably  provide  the  basic  elements  of  a  new  infrasonlc  monitoring 
network.  The  practical  implications  of  this  observation  are  worth  noting. 
Instead  of  developing  and  deploying  a  new  generation  of  microbarograph 
arrays,  it  may  be  possible  to  acquire  a  similar  capability  by  adding  mic-o- 
barographs  to  the  instrumentation  already  in  place  at  selected  sites  in  an 
existing  seismic  network.  The  cost  of  acquiring  an  infrasonlc  monitoring 
capability  in  this  manner  will  be  substantially  lower  than  the  deployment  and 
operation  of  an  independent  network  of  microbarograph  arrays.  The  basic 
objective  of  the  research  program  currently  underway  is  to  provide  the  tech¬ 
nical  basis  for  evaluating  seismograph-mlcrobarograph  options. 

The  basic  structure  of  the  research  program  may  be  deduced  from  the  statement 
of  work  shown  in  figure  2.  The  purpose  of  this  report  is  to  summarize 
progress  made  during  the  first  six  months  of  the  program. 


Task  1.  Acquire  a  combined  microbarometric  and  seismic  data  base 
pertaining  to  infrasonic  signals,  subsonic  noise,  and  seismic 
background  noise. 

Task  2.  Utilize  the  acquired  data  base  to  experimentally  assess  the  current 
capability  to  detect  and  identify  infrasonic  signals  with  combina¬ 
tions  of  the  outputs  of  a  microbarograph  and  a  three-component 
inertial  seismograph  system. 

Task  3.  Investigate  the  feasibility  of  acquiring  multi-component  earth 
strain  data  in  addition  to  microbarometric  and  inertial  seismic 
data  to  assist  in  the  detection  and  identification  of  infrasonic 
signal.  Also  investigate  modifications  to  existing  sensor  systems 
and  data  processing  scheme  which  would  contribute  to  the 
detection  and  identification  of  both  infrasonic  as  well  as  seismic 
signals. 


FIGURE  2.  SUMMARY  STATEMENT  OF  WORK 
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2.0  DATA  ACQUISITION 

The  primary  objectives  of  the  first  phase  of  the  research  program  was  to 
establish  a  temporary  observatory  for  the  simultaneous  acquisition  of  data 
pertaining  to  earth  movements  and  atmospheric  pressure  oscillations.  This 
objective  has  been  accomplished,  and  data  are  now  being  routinely  recorded 
from  a  five-element  microbarograph  array  and  a  three-component  KS-36000  long- 
period  seismograph  system  located  at  a  depth  of  150  meters  at  a  site  near 
McKinney,  Texas.  The  current  configuration  of  the  microbarograph  array  and 
its  location  with  respect  to  the  seismograph  system  is  shown  in  figure  3. 

The  microbarographs  are  of  the  M-4  generation  and  were  loaned  to  the  program 
by  S.M.U.  A  50-foot  length  of  garden  hose  with  0.75  inch  inside  diameter  Is 
attached  to  the  inlet  port  of  each  microbarograph.  Hypodermic  needles  are 
inserted  at  intervals  of  approximately  1  meter  to  provide  some  suppression  of 
very  short  wavelength  atmospheric  turbulence.  An  anemometer  is  also 
installed  near  the  borehole  containing  the  seismograph  system  to  provide  con¬ 
tinuous  information  on  the  local  surface  wind  speeds. 

The  seismograph  is  a  KS-36000  three-component  system  and  is  located  at 
33°14,56"N  latitude,  96°39,07”W  longitude  at  a  depth  of  152  meters.  Data 
from  all  sensing  units  are  transmitted  via  cable  to  a  central  data  collection 
facility  where  they  are  conditioned,  then  recorded  on  35mm  film  for  moni¬ 
toring  purposes.  In  addition,  the  data  are  sampled  at  0.5  second  intervals 
and  stored  on  digital  magnetic  tape  for  future  data  processing  and  analysis. 

A  block  diagram  of  the  data  acquisition  system  is  shown  in  figure  4.  The 
sensor  system  responses  are  shown  in  figure  5. 


SEISMIC  SYSTEM  / 

MCB  *  1 


FIGURE  3.  FIELD  CONFIGURATION  OF  THE  MICROBAROGRAPH  ARRAY  AND  THE 
DEEP  WELL  SEISMIC  INSTALLATION 
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FIGURE  4.  BLOCK  DIAGRAM  OF  SEISMIC  AND  RECORDING  SUBSYSTEMS 


FREQUENCY  (Hz) 


FIGURE  5.  RELATIVE  RESPONSES  OF  SEISMOGRAPH  CHANNELS  TO  DISPLACEMENT  AND 
MICROBAROGRAPH  CHANNELS  TO  PRESSURE 
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3.0  DATA  ANALYSIS 


3.1  PRESSURE-DISPLACEMENT  TRANSFER  FUNCTIONS 

The  passage  of  an  infrasonlc  wave  creates  quasi-static  earth  movements  which 
may  be  detected  at  the  output  of  a  seismograph  system.  The  amplitude  of 
these  earth  movements  may  be  related  to  the  amplitude  of  the  infrasonlc  wave 
through  a  transfer  function  whose  form  Is  dependent  upon  the  local  distribu¬ 
tion  of  elastic  constants,  the  speed  of  the  infrasonlc  wave  and  the  depth  of 
observation  (Sorrells,  1971).  By  approximating  the  earth  structure  at  the 
McKinney  observatory  as  a  layered  elastic  solid,  theoretical  transfer  func¬ 
tions  may  be  calculated  using  techniques  described  by  Sorrells  and  Goforth, 
1973.  Vertical  and  horizontal  transfer  functions  appropriate  for  an  obser¬ 
vational  depth  of  150  meters  at  McKinney,  Texas,  are  shown  In  figure  6.  The 
apparent  horizontal  displacement  transfer  function  is  the  sum  of  the  true 
horizontal  displacement  and  the  apparent  horizontal  displacements  caused  by 
earth  tilts.  It  is  Important  to  note  that  for  a  pressure  wave  propagating  at 
infrasonlc  speeds,  the  vertical  displacements  will  generally  be  greater  than 
the  apparent  horizontal  displacements  at  frequencies  higher  than  about  0.006 
Hz.  As  shown  in  figure  7,  this  crossover  point  shifts  rapidly  to  higher  fre¬ 
quencies  as  the  propagation  speed  of  the  pressure  disturbance  decreases. 

Now,  generally  speaking,  the  propagation  speed  of  wind-related  turbulence  is 
approximately  equal  to  the  mean  wind  speed  tdiich  rarely  exceeds  10  m/sec  for 
any  sustained  length  of  time.  The  comparisons  of  vertical  and  horizontal 
earth  movements  provide  the  basis  for  a  rapid  means  of  roughly  classifying 
pressure-related  earth  movements  into  potential  signal  and  noise  categories; 
i.e.,  if  the  observed  horizontal  component  of  earth  movement  is  greater  than 
the  observed  vertical  component  and  its  frequency  is  greater  than  0.006  Hz, 
then,  the  disturbances  can  be  safely  classified  as  noise.  If  these  criteria 
are  not  satisfied,  then  the  disturbance  has  a  high  probability  of  being 
infrasonlc  in  nature.  This  technique  is  now  being  Implemented  to  select  data 
records  for  further  processing  and  analysis. 
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FIGURE  6.  RESPONSE  OF  AN  ELASTIC  LAYERED  HALFSPACE  TO  ATMOSPHERIC 
LOADING  AT  INFRASONIC  VELOCITY 
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FIGURE  7.  RESPONSE  OF  AN  ELASTIC  LAYERED  HALFSPACE  TO  ATMOSPHERIC 
LOADING  AT  SUBSONIC  VELOCITY 
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3.2  INFRASONIC  SNR'S:  PRELIMINARY  RESULTS 


The  primary  objective  of  the  initial  phase  of  this  experiment  was  to 
establish  a  temporary  facility  for  the  simultaneous  acquisition  of  digital 
data  pertaining  to  local  earth  movements  and  atmospheric  pressure 
oscillation.  This  objective  has  been  achieved,  and  we  are  now  in  the  process 
of  experimentally  evaluating  the  combined  use  of  seismic  and  microbarometric 
data  to  detect  and  identify  lnfrasonlc  signals.  Preliminary  results  bearing 
on  this  subject  are  briefly  summarized  in  the  following  paragraphs. 

It  was  shown  in  an  earlier  section  that  it  was  theoretically  possible  to 
obtain  a  significant  improvement  in  the  lnfrasonlc  signal -to-wind  noise  ratio 
by  transforming  the  detection  problem  into  the  seismic  regime.  In 
particular,  it  was  shown  that  because  of  the  passive  wavelength  filter  pro¬ 
perties  of  the  earth  with  respect  to  atmospheric  pressure  variations,  the 
lnfrasonlc  signal-to-wlnd  noise  ratio  must  always  be  greater  at  the  output  of 
a  vertical  seismograph  than  at  the  output  of  a  microbarograph.  In  practice, 
however,  the  actual  gain  in  SNR  obtained  by  utilizing  a  vertical  seismograph 
system  will  be  limited  by  the  existence  of  earth  noise  that  is  unrelated  to 
atmospheric  pressure  variations.  At  the  present  time,  there  is  little 
experimental  evidence  bearing  upon  this  subject.  Thus,  a  first  objective  of 
the  data  analysis  is  to  quantitatively  assess  the  nature  of  this  limitation. 
The  approach  currently  being  followed  is  outlined  below. 

Let  m  (t)  and  uz  (t)  denote  the  outputs  of  a  microbarograph  and  a  vertical 
seismograph,  respectively.  Then,  when  the  wind  is  blowing  and  an  lnfrasonlc 
signal  is  present: 

m(t)  -  rm  *  (p8  +  pn)  (1) 

Where  rm  is  the  response  of  the  microbarograph,  p8  is  the  lnfrasonlc  pressure 
signal,  pn  is  the  turbulent  pressure  noise  created  by  the  wind  and  (*)  deno¬ 
tes  convolution.  Similarly, 


-  rz  *(gz8  *  p8  +  gzn  *  Pn  +  nz> 


(2) 


Where  gz8  and  gzn  are  the  Greens  functions  for  the  vertical  component  of 
displacement  caused  by  lnfrasonlc  signal  and  wind  noise,  respectively.  For  a 
vertical  seismograph  located  at  a  depth  of  150  meters,  it  can  be  shown  that 


«8zs  *  Pn>2>  «  «nz>2>  (3) 

so  that  for  our  purposes,  equation  2  may  be  approximated  as: 

uz(t)  =  rz*(gz8  *  p8  +  nz)  (4) 

Now,  the  corresponding  power  spectra  may  be  written  as: 

M(w)  -  RapCPs  +  Pn)  (5) 

Uz(w)  =  Rz  |2  (  |Gzs  |2  ps  +  Nz>  <6> 


Where  M,  Uz,  Pg,  Pn  and  Ns  are  the  power  spectral  density  functions  associated 
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with  m,  Ug,  p8,  pn  and  nj,  and  Rq,  and  Gzs  are  the  complex  terms  for  func¬ 
tions  associated  with  rm,  rz  and  gz8.  Rewriting  equation  6  as: 


r*.  wr.  i2  =  Ps  + 


~  p  +  p ' 

ra  r  n 


and  rewriting  equation  5  as: 
M 


“  *s  + 


it  can  be  seen  that  the  term 
N0 


P'n  - 


8 


(7) 


(8) 


(9) 


can  be  thought  of  as  an  equivalent  pressure  noise  power  spectral  density 
estimate  for  the  earth  motion  field. 


Let  (SNR)Z  and  (SNR)m  be  the  signal-to-nolse  ratios  for  the  earth  motion  and 
pressure  fields,  then 


(SNR)Z  -  -£=■ 

and 

(SNR)m  - 


The  ratio 


I 


zm 


(SNR)Z 

(SNR)ra 


(10) 


(ID 


(12) 


is  thus  a  measure  of  the  actual  Improvement  In  infrasonic  signal -to-noise 
ratios  that  can  be  obtained  through  the  use  of  a  vertical  seismograph. 
Notice  that 


Thus,  if  we  confine  our  attention  to  Intervals  triten  Pg  ~  0  and,  if  we  know 
the  moduli  of  the  transfer  functions  Rg,,  Rj  and  Gsz,  we  can  obtain  an  esti¬ 
mate  of  Izn  from  the  outputs  of  a  microbarograph  and  a  vertical  seismograph. 
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We  are  currently  making  estimates  of  Izm  during  various  states  of  atmospheric 
turbulence  given  the  assumption  that  the  moduli  of  the  transfer  functions 
shown  in  figures  5  and  6  accurately  depict  the  instrumental  and  earth 
response  functions.  An  important  preliminary  result  is  illustrated  in 
figures  8a  and  8b.  Here,  we  compare  estimates  of  Pn(f)  and  P'n(f)  in  a  band 
extending  from  approximately  0.002  to  0.5  Hz  (figure  8a)  and  display  the 
corresponding  estimate  of  Izm(f)  (figure  8b).  These  results  were  obtained 
from  spectral  estimates  of  outputs  of  the  vertical  seismograph  and  a  nearby 
surface  microbarograph  during  an  interval  of  gusty  winds.  The  spectral  esti¬ 
mates  were  made  by  applying  the  Welch  method  of  averaging  periodograms  to 
time  series  consisting  of  10  blocks  of  1024  points.  The  important  point 
illustrated  by  these  results  is  that  the  estimated  value  of  P'n(f)  is  signi¬ 
ficantly  lower  than  the  estimates  of  Pn(f)  over  a  frequency  range  extending 
from  approximately  0.005  to  0.05  Hz.  This  means,  in  effect,  that  within  this 
frequency  band,  the  infrasonlc  SNR  will  be  greater  at  the  output  of  a  ver¬ 
tical  seismograph  than  at  the  output  of  a  microbarograph  during  this 
interval.  The  estimate  of  Izm(f),  which  is  shown  in  figure  8b,  indicates 
that  the  increase  in  infrasonlc  SNR's  can  be  as  much  as  15  dB  in  a  narrow 
frequency  band  near  0.03  Hz.  These  results  are  encouraging  in  that  they  are 
qualitatively  consistent  with  earlier  observations  made  at  Grand  Saline, 
Texas,  (Sorrells  et  al,  1971)  which  Indicated  that  over  a  similar  frequency 
range,  the  infrasonlc  SNR's  for  a  presumed  atmospheric  nuclear  explosion  were 
higher  at  the  output  of  a  vertical  seismograph  located  at  a  depth  of  approxi¬ 
mately  180  meters  than  at  the  output  of  a  microbarograph  located  at  the 
surface . 

Similar  calculations  have  been  made  for  data  collected  during  calm  Intervals. 
These  results  are  shown  in  figures  9a  and  9b.  Observe  that,  under  these 
circumstances,  Pn(f)  is  less  than  P'n(f)  throughout  the  entire  bandwidth  of 
interest.  These  results  simply  serve  to  illustrate  the  obvious  point  that  in 
the  absence  of  wind-generated  atmospheric  pressure  changes,  a  microbarograph 
will  be  superior  to  a  vertical  seismograph  for  the  detection  of  infrasonlc 
waves . 
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(a) 


(b) 
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FIGURE  9.  (a)  COMPARISON  OF  OBSERVED  (SOLID  LINE)  AND  SEISMIC  EQUIVALENT  (DASHED 
LINE)  PRESSURE  NOISE  SPECTRAL  ESTIMATES,  (b)  ESTIMATED  INFRASONIC  SNR 
ENHANCEMENT.  DATA  SHOWN  WERE  COLLECTED  DURING  A  CALM  INTERVAL. 
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4.0  FUTURE  RESEARCH  DIRECTIONS 


The  preliminary  results  presented  in  the  previous  section  increase  our  con¬ 
fidence  that,  during  intervals  of  atmospheric  turbulence,  the  infrasonic  SNR 
will  be  greater  at  the  output  of  a  vertical  seismograph  than  at  the  output  of 
a  microbarograph  over  a  significant  fraction  of  the  frequency  range  of 
interest  for  atmospheric  nuclear  test  detection.  While  these  results  show 
promise,  more  work  needs  to  be  done  not  only  to  document  the  gain  in  the 
infrasonic  SNR  both  as  a  function  of  frequency  and  the  intensity  of 
atmospheric  turbulence,  but  also  to  refine  our  estimates.  In  particular, 
experimental  confirmation  of  predicted  pressure-displacement  transfer  func¬ 
tion  is  needed  since  it  plays  a  significant  role  in  estimation  of  the  SNR 
gain. 

It  is  also  important,  for  the  purposes  of  this  program,  to  renew  investiga¬ 
tions  into  the  relationship  between  the  yield  of  an  atmospheric  nuclear 
explosion  and  the  spectrum  of  the  infrasonic  signal  vAiich  it  generates. 

While  it  is  generally  recognized  that  the  infrasonic  spectrum  shifts  toward 
higher  frequencies  as  the  yield  decreases,  the  precise  nature  of  the  rela¬ 
tionship  is  poorly  known,  particularly  in  the  kiloton  yield  range.  A 
knowledge  of  this  relationship  is  important  in  that  it  would  allow  us  to 
translate  the  observed  frequency  limits  currently  being  determined  to  a  yield 
range.  Be  knowing  the  yield  range  over  which  enhanced  infrasonic  SNR's  are 
possible,  the  technical  feasibility  of  developing  new  technology  to  extend 
the  currently  observed  frequency  range  can  be  realistically  assessed. 

Finally,  the  program  should  provide  a  reliable  estimate  of  the  infrasonic 
signal  detection  threshold  for  a  combination  of  sensors  consisting  of  a 
three-component  seismograph  and  a  microbarograph.  It  is  probable  that  this 
threshold  can  be  determined  directly  from  the  experimental  data  base  for  fre¬ 
quencies  lower  than  about  0.02  Hz  utilizing  naturally  occurring  infrasonic 
signals.  At  higher  frequencies,  an  indirect  method  of  estimation  may  be 
needed . 
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